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Abstract 

The  paper  presents  results  of  computer  experiments  per¬ 
formed  for  research  to  enhance  operational  parameters  of  multi¬ 
cavity  magnetrons  through  using  the  metal-thin-wire  (MTW) 
and  the  metal-split-ring  (MSR)  metamaterial  structures  as  a  cold 
(non-thermionic)  cathode.  It  is  demonstrated  by  particle-in-cell 
(PIC)  simulations  of  a  generic  10-vane  non-relativistic  ultra-high 
frequency  (UHF)  magnetron  that  the  output  microwave  power 
and  the  start-up  time  of  this  magnetron  are  improved  by  replac¬ 
ing  either  a  solid  or  a  transparent  cathode  by  a  bulk  MTW  me¬ 
tamaterial  structure  formed  by  an  array  of  individual  wires  di¬ 
rected  parallel  to  the  axis  of  the  magnetron  (rodded  cathode). 

I.  INTRODUCTION 

The  artificial  MTW  and  the  MSR  structures  possess 
metamaterial  properties,  i.e.,  negative-e/positive-/i  and 
positive-£/negative-/i,  respectively,  when  the  incident 
time-varying  electric  field  Ex  oscillates  along  the  wires  of 
the  MTW  structure,  or  when  the  incident  time-varying 
magnetic  Hx  field  oscillates  perpendicularly  to  the  rings  of 
the  MSR  structure  [1],  [2].  In  both  cases,  the  incident  Ex 
and  ELX  fields  induce  electric  currents  flowing  along  the 
individual  metal  wires  (rings).  The  oscillations  of  the  in¬ 
duced  currents  generate  the  artificial  electric  or  magnetic 
dipole  moments  within  the  bulk  of  a  metamaterial  struc¬ 
ture,  which  either  constructively  or  destructively  interact 
with  the  incident  fields.  Both  MTW  and  MSR  structures 
may  be  considered  as  quasi-metamaterials  (transparent) 
structures  when  the  incident  time -varying  electric  field  Ex 
is  orthogonal  to  the  individual  wires  and  freely  penetrates 
inside  of  those  structures. 

The  multi-cavity  magnetron  operates  in  one  of  the 
transverse-electrical-//^  ( TE-like )  cavity  modes  of  the 
magnetron  resonant  cavity  (magnetron  diode)  formed  by 
the  cathode  emitting  electrons  and  the  anode  consisting  of 
a  number  of  resonators  [3],  [4].  The  resonant  cavity  slows 
down  the  induced  electromagnetic  wave  azimuthally  trav¬ 
eling  between  the  cathode  and  the  anode.  The  spatial  dis¬ 
tribution  of  the  time-varying  electromagnetic  field 
( ExxH i),  induced  within  the  resonant  cavity  and  asso¬ 
ciated  with  the  traveling  electromagnetic  wave,  is  charac¬ 
terized  by  the  circumferential  Ex§  and  the  radial  EXp  com¬ 
ponents  of  the  electric  field  and  the  axial  HXz  component 
of  the  magnetic  field.  The  number  of  Ex§  phase  variations 
in  the  azimuthal  direction  is  determined  by  the  operating 
frequency  of  the  magnetron  (mode  of  the  magnetron  op¬ 
eration)  and  the  anode  geometry  (number  of  cavities).  The 
number  of  EXp  and  HXz  phase  variation  along  the  radius 
and  the  axis  of  the  resonant  system  are  usually  minimized 
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to  one.  This  allows  the  magnetron  to  operate  in  the  lower- 
frequency  TE-like  mode  (at  the  given  number  of  Ex§  cir¬ 
cumferential  variations)  and  prevents  axial  and  radial 
mode  competition.  The  external  electric  field  E0  in  the 
magnetron  diode  is  directed  radially  from  the  anode  to  the 
cathode.  The  external  magnetic  field  H0  is  directed  axial¬ 
ly.  The  crossed  E0xH0  fields  force  electrons  emitted  from 
the  cathode  to  drift  in  the  azimuthal  direction.  The  multi - 
cavity  magnetron  starts  to  generate  (induce)  ExxHx  fields 
and  produce  microwaves  (output  power)  when  the  drift 
velocity  of  the  emitted  electron  is  approximately  equal  to 
the  phase  velocity  of  Ex§. 

For  the  multi-cavity  magnetron  operating  in  the  TE-like 
mode,  the  bulk  metamaterial  cathode  may  be  designed  in 
the  form  of  a  single/double/multi-helix  structure  (Fig.  la), 
the  individual  elements  (turns)  of  which  are  oriented  pa¬ 
rallel  to  Eh  and  perpendicular  to  HXz.  The  proper  choos¬ 
ing  of  the  helix  pitch  may  allow  implementing  the  nega- 
tive-e/negative-/,*  properties  of  the  metamaterial  cathode 
within  the  desired  microwave  frequency  range  to  enhance 
and  improve  the  magnetron  operational  parameters.  The 
metamaterial  cathode  may  be  also  designed  in  the  form  of 
the  rodded  medium  whose  individual  elements  (rods)  are 
placed  either  radially  (Fig.  lb)  (parallel  to  EXp  and  perpen¬ 
dicular  to  both  Ex§  and  HXz)  of  longitudinally  (Fig.  lc) 
(perpendicular  to  both  EXp  and  Ex§  and  parallel  to  HXz).  In 
the  latter  case  (Fig.  lc),  the  rodded  medium  is  essentially  a 
quasi-metamaterial  structure. 


Fig.  1.  The  different  cathodes:  a)  multi-helix  structure  with  met¬ 
al  wires  parallel  to  Ex§  and  the  metal  rings  perpendicular  to  HXz, , 
b)  multi-point  structure  with  metal  wires  parallel  to  EXp  (multi¬ 
point  cathode),  c)  multi-rod  structure  with  metal  wires  parallel 
to  HXz  (rodded  cathode),  d)  standard  helical  cathode  with  non¬ 
emitting  central  rod. 
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Results  of  the  PIC  simulations  of  a  generic  10-vane 
non-relativistic  continuous -wave  (CW)  UHF  magnetron 
with  a  quasi-metamaterial  multi-rod  or  rodded  cathode 
(Fig.  lc)  are  described  below.  Most  important  geometrical 
dimensions  of  the  magnetron  are  similar  to  the  commer¬ 
cially  available  CWM-75/100L  and  L-4031  tubes  of  the 
L-3  Communications,  Electron  Devices  [5]. 

II.  SIMULATION  GEOMETRY  AND  PROCEDURE 

The  simulations  are  performed  using  ICEPIC  code  de¬ 
veloped  and  maintained  by  the  AFRL/RDHE  [6]-[8]. 
ICEPIC  is  a  hilly  relativistic,  three-dimensional  (3D), 
variable-mesh  PIC  code  capable  of  simulating  the  inter¬ 
disciplinary  physics  of  charge -particle  beams,  high-power 
microwaves  (HPM),  and  plasmas  in  complex  geometries 
of  modern  microwave  vacuum  electronic  devices 
(MVEDs)  [8].  The  code  incorporates  numerous  advances 
in  parallel  computing  and  performs  automatic  partitioning 
using  a  single  domain  decomposition  algorithm,  asyn¬ 
chronous  message  passing,  and  dynamic  load  balancing 
[7].  ICEPIC  simulates  from  the  first  principles  (Maxwell’s 
equations  and  Lorentz’s  forces  law)  the  relativistic  elec¬ 
trodynamics  and  charged-particle  dynamics  within  the 
resonant  structures  of  different  MVEDs  [6]. 

The  simulation  model  of  the  magnetron  is  built  within 
the  Cartesian  x-y-z  geometry  using  simulation  cells  with 
grid  size  of  0.5x0. 5x0. 5  mm3  [9].  It  consists  of:  (i)  the 
cathode  block  with  the  input  port,  (ii)  the  10-vane  anode 
block  surrounding  the  emitting  part  of  the  cathode  block, 
(iii)  the  “double-ring-strapping”  system  coupling  the  al¬ 
ternative  resonators  of  the  anode  block,  and  (iv)  three 
output  electrodes  connecting  three  vanes  of  the  anode 
block  with  the  appropriate  output  ports. 

The  anode  block  consists  of  ten  vanes,  7V=10,  forming  a 
cylindrical  slow- wave  structure  (SWS).  There  are  two 
grooves,  one  on  each  end  of  the  anode  block.  The  double¬ 
ring-strapping  system  is  located  within  the  grooves  on 
each  end  of  the  anode  block.  Each  strap  electrically  con¬ 
nects  alternate  vanes  and  passes  over  adjacent  vanes  of 
the  anode  block.  Three  cylindrical  electrodes  connect  the 
first,  the  fourth,  and  the  sixth  vanes  of  the  anode  block 
with  the  output  ports.  The  cathode  block  consists  of  the 
non-emitting  central  electrode  connecting  the  input  port 
with  the  emitting  part  of  the  cathode,  two  non-emitting 
end  caps  preventing  emitted  electrons  from  leaking  out 
the  interaction  space,  and  the  cathode  itself  (Fig.  2). 

The  simulations  of  the  magnetron  are  performed  at  in¬ 
put  voltage  Fin=45  kV  and  magnetic  flux  density  B^=0A9 
T.  The  rise  time  of  Vm  is  10  ns  [10]  at  simulation  time  of 
2000  ns.  This  allows  one  to  consider  Vm  as  the  “steady- 
state”  input  parameter.  B 0  is  axially  directed,  uniform,  and 
constant  within  the  simulation  volume  and  during  the  si¬ 
mulation  time. 

In  the  course  of  the  simulations,  Vin  is  monitored  (as 
well  as  the  input  current  7in  and  the  input  power  Pin)  at  the 
input  port.  The  output  power  Poui  is  monitored  (as  well  as 
the  output  rf  voltage  Vout)  at  the  output  ports.  The  resona¬ 


tor  voltage  Vres  is  monitored  inside  the  first  resonator  of 
the  anode  block.  The  anode  current  /a  is  monitored  at  the 
plane  located  on  a  small  distance  (one  or  two  simulation 
cells)  from  the  anode  surface.  All  monitored  parameters 
are  smoothed  over  100,000  consecutive  time  steps  and 
averaged  over  their  steady-state  phases  (when  the  moni¬ 
tored  parameter  is  not  changed  in  time)  in  order  to  get  the 
steady-state  characteristics  of  the  magnetron  operation. 

The  entire  simulation  procedure  develops  as  follows. 
Once  VXn  is  applied  to  the  input  port,  the  space -charge 
limited  (SCL)  algorithm  of  the  ICEPIC  code  [11]  begins 
to  produce  electrons  at  a  small  distance  from  the  cathode 
(one  or  two  grid  cells).  These  originated  electrons  are 
extracted  from  the  cathode  and  accelerated  by  E0 ,  def¬ 
lected  by  H0=juoB0  (juq  is  permeability  of  vacuum),  and 
move  within  the  interaction  space  in  accordance  with  the 
given  distributions  of  E0xH0  and  the  electron  space- 
charge  qe  accumulated  within  the  interaction  space. 

The  accelerated  electrons  give  rise  to  Im  flowing  along 
the  central  electrode  connecting  the  input  port  with  the 
cathode.  After  some  amount  of  time,  determined  by  the 
loaded  Q  (quality  factor  [4])  of  the  magnetron  resonant 
cavity,  the  magnetron  spokes  begin  to  grow  and  the 
magnetron  starts  to  operate.  After  formation  of  the 
magnetron  spokes,  both  /a  and  Pout  increase  and  go  into 
their  steady-state  phases.  The  electrons  begin  to  move 
within  the  interaction  space  in  accordance  with  the  given 
distributions  of  E0xH0 ,  ExxHx ,  and  qe. 

Results  of  the  simulations  of  the  magnetron  having  the 
rodded  cathode  (Fig.  lc,  Fig.  2c)  are  compared  with  the 
results  obtained  during  PIC  simulations  of  the  same 
magnetron  having  the  solid  cylindrical  (Fig.  2a)  and  the 
transparent  [12],  [13]  (Fig.  2b)  cathodes.  All  cathodes  of 
the  simulation  model  of  the  magnetron  (Fig.  2)  are  sup¬ 
posed  to  operate  in  the  explosive  electron  emission  (EEE) 
mode  [14],  which  automatically  assumes  the  space-charge 
limited  (SCL)  mechanism  of  the  electron  extraction  from 
the  cathode  into  the  interaction  space.  Note,  however,  that 
the  standard  high-power  industrial-heating  UHF  CWM- 
75/1 00L  and  L-4031  MHz  magnetrons  have  the  helical 
thermionic  cathode  (with  the  central  back-current  elec¬ 
trode)  (Fig.  Id)  working  in  the  thermal  electron-emission 
(TEE)  mode  [5]. 

III.  SIMULATION  RESULTS 

The  particle  plot  distributions  within  the  interaction 
space  of  the  magnetron  having  three  different  cathodes 
(Fig.  2)  (solid,  transparent,  and  rodded)  are  shown  in  Fig.  3 
(-10  ns  of  simulation  time),  Fig.  4  (-100  ns),  and  Fig.  5 
(-1000  ns).  It  can  be  seen  in  Fig.  -  ig.  5  that  the  magne¬ 
tron  begins  to  oscillate  in  the  7i-mode  with  either  cathode. 
The  7r-mode  is  recognized  by  five  magnetron  spokes  cor¬ 
responding  to  the  TE5X  mode  of  ExxHx  distribution  within 
the  interaction  space. 

There  are,  however,  very  significant  differences  in  the 
amount  of  the  anode  current  /a  and  the  output  power  Pout 
when  the  magnetron  operates  with  the  different  cathodes. 
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Fig.  2.  Cathodes  used  in  simulations:  a)  solid  cathode,  b)  transparent  cathode,  c)  rodded  (quasi-metamaterial)  cathode  (Fig.  lc). 


a  b  c 

Fig.  4.  Particle  plots  at  100  ns  of  simulation  time:  a)  solid  cathode,  b)  transparent  cathode,  c)  rodded  (quasi-metamaterial)  cathode. 


Fig.  5.  Particle  plots  at  1000  ns  of  simulation  time:  a)  solid  cathode,  b)  transparent  cathode,  c)  rodded  (quasi-metamaterial)  cathode. 
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a  b  c 

Fig.  6.  Input  Voltage,  Vm:  a)  solid  cathode,  b)  transparent  cathode,  c)  rodded  (quasi -metamaterial)  cathode. 


Fig.  8.  Output  Power,  Pout :  a)  solid  cathode,  b)  transparent  cathode,  c)  rodded  (quasi -metamaterial)  cathode. 
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Fig.  10.  Spectrum  of  Output  Voltage,  Vout :  a)  solid  cathode,  b)  transparent  cathode,  c)  rodded  (quasi -metamaterial)  cathode. 
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The  differences  can  be  seen  in  Fig.  7  that  the  steady- 
state  /a  increases  from  -7.7  A  in  the  magnetron  having 
solid  cathode,  to  -16.7  A  in  the  magnetron  having  trans¬ 
parent  cathode,  and  up  to  -23 . 1  A  in  the  magnetron  hav¬ 
ing  rodded  cathode.  In  the  exact  accordance  with  4  in¬ 
crease  [15],  [16],  the  steady-state  Pout  increases  from  -311 
kW  in  the  magnetron  having  solid  cathode,  to  -62 1  kW  in 
the  magnetron  having  transparent  cathode,  and  up  to  -820 
kW  in  the  magnetron  having  rodded  cathode  (Fig.  8). 

Detailed  analysis  of  the  obtained  results  show  that  both 
4  (Fig.  7)  and  Pout  (Fig.  8)  increase  in  direct  proportion  to 
the  square  of  the  resonator  voltage  Vres  (Fig.  9),  which  is 
directly  proportional  to  the  induced  time -varying  electric 
field  amplitude  Ex.  The  frequency  of  the  magnetron  op¬ 
eration  fn  is  determined  by  the  frequency  of  Ex  oscilla¬ 
tions.  Frequency  fn  remains  the  same  with  all  cathodes 
(Fig.  10). 

The  numerical  relations  between  4,  Pout  and  Vres  dem¬ 
onstrate  that  both  4  and  Pout  increase  in  direct  proportion 
to  the  square  of  Fres.  This  allows  one  to  suppose  that  the 
direct  cause  of  both  4  and  Pout  increase  is  nothing  but  the 
E i  increase  within  the  interaction  space.  The  smoothed 
traces  of  4  and  Pout  are  shown  in  Fig.  11  and  Fig.  12,  cor¬ 
respondent^.  The  envelope  of  Vres  is  shown  in  Fig.  13.  All 
traces  demonstrate  how  the  cathode  construction  (Fig.  2) 
affects  the  magnetron  operation. 
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Fig.  11.  Smoothed  traces  of  the  anode  current  Ia. 
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Fig.  12.  Smoothed  traces  of  the  output  power  Pout- 


Fig.  13.  Envelope  of  the  resonator  voltage  Fres. 

IV.  CONCLUSION 

PIC  simulations  of  a  generic  10-vane  UHF  magnetron 
with  the  solid,  the  transparent,  and  the  quasi -metamaterial 
(rodded)  cathodes  (Fig.  2)  are  performed.  Results  of  the 
simulation  showed  that  the  use  of  the  rodded  cathode  al¬ 
lows  significantly  increasing  the  output  microwave  power 
(Fig.  12)  by  the  appropriate  increase  of  the  anode  current 
(Fig.  11)  of  the  magnetron.  The  anode  current  increase  is 
caused  by  appropriate  increase  of  the  time -varying  in¬ 
duced  electric  field  Ex  (square  of  the  resonator  voltage 
(Fig.  13))  within  the  interaction  space.  The  latter  effect  is 
the  direct  result  of  the  metamaterial  structure  (rodded  ca¬ 
thode)  use  within  the  magnetron  resonant  cavity. 

V.  REFERENCES 

[1]  L.  Solimar  and  E.  Shamonina,  Waves  in  Metamaterial , 
Oxford  University  Press,  2009. 

[2]  C.  Caloz  and  T.  Itoh,  Electromagnetic  Metamaterials: 
Transmission  Line  Theory  and  Microwave  Applications , 
John  Wiley  &  Sons,  Inc.,  2006. 

[3]  R.M.  Gilgenbach,  Y.Y.  Lau,  H.  McDowell,  K.L. 

Cartwright,  and  T.A.  Spencer,  “Crossed-field  devices,”  in 
Modern  Microwave  and  Millimeter-  Wave  Power  Electron¬ 
ics  Edited  by  R.J.  Barker,  J.H.  Booske,  N.  C.  Luhmann,  Jr, 
and  G.S.  Nusinovich ,  IEEE  Press  and  John  Wiley  &  Sons, 
Inc.,  2005,  pp.289-342. 

[4]  J.  Benford,  J.A.  Swegle,  and  E.  Schamiloglu,  “ Relativistic 
magnetrons  and  MILOs ,”  in  High-Power  Microwaves, 
Second  Edition ,  CRC  Press  and  Taylor  &  Francis  Group, 
2007,  pp.259-320. 

[5]  Wynn,  A.P.,  D.E.  Blank,  P.S.  Campbell,  R.R.  Lentz,  W.T. 
Main,  S.G.  Tantawi,  K.G.  Kato,  H.K.  Beutel,  K.W.  Brown, 
D.D.  Crouch,  G.K.  Jones,  and  R.B.  McDonald,  “Develop¬ 
ment  of  a  300  kW  CW  L-band  industrial  heating  magne¬ 
tron,”  in  Fifth  IEEE  International  Vacuum  Electronics 
Conference  IVEC  2004,  April  27-29,  2004,  pp.  164- 165. 

[6]  G.E.  Sasser,  L.  Bowers,  S.  Colella,  D.  Lileikis,  J.  Lugins- 
land,  D.  McGrath,  L.  Merkle,  R.E.  Peterlcin,  and  J.J.  Watr- 
ous,  “Virtual  prototyping  of  microwave  devices  using 
MIID.  PIC,  and  CEM  codes,”  Proceedings  of  the  29th  AIAA 
Plasmadynamics  and  Lasers  Conference,  Albuquerque, 

NM,  USA,  15-18  June  1998 ,  publication  #  AIAA-98-2926, 
AIAA  Press,  1998. 


Approved  for  public  release;  distribution  is  unlimited. 


[7]  K.L.  Cartwright,  A.G.  Greenwood,  P.J.  Mardahl,  T.  Mur¬ 
phy,  and  M.D.  Haworth,  “Virtual  prototyping  of  radio  fre¬ 
quency  weapons,”  Journal  of  Directed  Energy ,  vol.  1,  no. 

1,  Fall  2003,  pp.  16-34. 

[8]  D.  Shiftier,  T.K.  Staton,  T.W.  Hussey,  O.  Zhou,  and  P. 
Mardahl,  “High-power  microwave  sources,”  in  Modern 
Microwave  and  Millimeter-Wave  Power  Electronics ,  edited 
by  R.J.  Barker,  J.H.  Booske,  N.C.  Luhmann  Jr.,  and  G.S. 
Nusinivich,  John  Wiley  &  Sons,  2005,  pp.691-730. 

[9]  A.D.  Greenwood,  “An  ICEPIC  convergence  study  using  a 
relativistic  magnetron,”  Air  Force  Research  Laboratory 
Technical  Memorandum ,  AFRL-DE-PS-TM-2005-1005, 
2005. 

[10]  D.  Price,  E.  Goldman,  D.  Nett,  S.  Putnam,  D.  Weidenhei- 
mer,  C.  Bloemker,  R.  White,  and  T.  Wynn,  “DEW  technol¬ 
ogy  at  the  Titan  Corporation,”  Proceedings  of  the  1st  AIA A 
International  Energy  Conversion  Engineering  Conference, 
Portsmouth,,  Virginia,  USA,  17-21  August  2003,  publica¬ 
tion  #  AIAA-2003-5913,  AIAA  Press,  2003. 

[11]  J.J.  Watrous,  J.W.  Luginsland,  and  G.E.  Sasser  III,  “An 
improved  space-charge-limited  emission  algorithm  for  use 
in  particle-in-cell  codes,”  Physics  of  Plasmas,  vol.  8,  no.  1, 
pp.  289-296,  January  2001. 

[12]  M.  Fulcs  and  E.  Schamiloglu,  “Rapid  start  of  oscillations  in 
a  magnetron  with  a  “transparent”  cathode,”  Physical  Re¬ 
view  Letters,  vol.  95,  no.  20,  1 1  November  2005,  pp. 
205101(1-4). 

[13]  M.C.  Jones,  V.B.  Neculaes,  Y.Y.  Lau,  R.M.  Gilgenbach, 
W.M.  White,  B.W.  Hoff,  and  N.M.  Jordan,  “Magnetron 
priming  by  multiple  cathodes,”  Applied  Physics  Letters , 
vol.  87,  2005,  pp.  081501(1-3). 

[14]  G.A.  Mesyats,  Explosive  Electron  Emission ,  URO  Press, 
1998. 

[15]  R.B.  Miller,  “The  relativistic  microwave  magnetron,”  in  An 
Introduction  to  the  Physics  of  Intense  Charged  Particle 
Beams ,  Plenum  Press,  1982,  pp.214-236. 

[16]  A.D.  Andreev,  K.J.  Hendricks,  M.I.  Fuks  and  E.  Schami¬ 
loglu,  “Elemental  theory  of  a  relativistic  magnetron  opera¬ 
tion:  Anode  current,  Journal  of  Directed  Energy ,  vol.  3,  no. 
4,  Summer  of  2010,  pp.349-383. 


Approved  for  public  release;  distribution  is  unlimited. 


DISTRIBUTION  LIST 


DTIC/OCP 

8725  John  J.  Kingman  Rd,  Suite  0944 
Ft  Belvoir,  YA  22060-6218  1  cy 

AFRL/RVIL 

Kirtland  AFB,  NM  87117-5776  2  cys 

Official  Record  Copy 

AFRL/RDHP/Kyle  Hendricks  1  cy 


